getwd()
setwd ("C:/Users/jess eduardo ortega/Documents/CURSO-R")
Eco <- read.table ("Nueva.csv", header=T, sep=";")
head (Eco)
str (Eco)
Ecofisiologia <- function(Eco)
  {
  if(class(Eco)!="data.frame") stop("\t", "Erro: O objeto x no  da classe numeric!\n") # Os dados devem ser numricos
  if(any(is.na(Eco)==TRUE)) stop("\t", "Erro: Os Objetos  de Temperaturas presenta NAs! Elimine os NAs correspondente a cada objeto y\n") # advertncia para eliminar NAs
  
########################## Cargando a funo Ecofisiologia ######################################

  x11()# Abre uma janela para grficos
  par(mfrow=c(2,2))# divide a janela em 2 fila e 2 colunas
  par(bty="l.5", tcl=-0.3)# me permite fazer o fechamento da caixa do grfico, muda a disposio da vinheta nos eixos 
  hist(Eco$Temperaturas[Eco$Variables=="Modelo"], main="Modelo", xlab="Temperatura  C", ylab="Frequncias", prob= T)# Histograma
  lines(density(Eco$Temperaturas [Eco$Variables=="Modelo"]),col="red")# mostra a linha de distribuio normal
  hist(Eco$Temperaturas[Eco$Variables=="Sp1"], main="Sp1", xlab="Temperatura  C", ylab="Frequncias", prob= T)# Histograma
  lines(density(Eco$Temperaturas [Eco$Variables=="Sp1"]),col="red")# mostra a linha de distribuio normal
  hist(Eco$Temperaturas[Eco$Variables=="Sp2"], main="Sp2", xlab="Temperatura  C", ylab="Frequncias", prob= T)# Histograma
  lines(density(Eco$Temperaturas [Eco$Variables=="Sp2"]),col="red")# mostra a linha de distribuio normal
  hist(Eco$Temperaturas[Eco$Variables=="Sp3"], main="Sp3", xlab="Temperatura  C", ylab="Frequncias", prob= T)# Histograma
  lines(density(Eco$Temperaturas [Eco$Variables=="Sp3"]),col="red")# mostra a linha de distribuio normal
  x11()# Abre uma janela para grficos
  par(mfrow=c(2,2))# divide a janela em 2 fila e 2 colunas
  par(bty="l.5", tcl=-0.3)# me permite fazer o fechamento da caixa do grfico, muda a disposio da vinheta nos eixos
  qqnorm(Eco$Temperaturas [Eco$Variables=="Modelo"],pch=16, col="black",main="Ajuste da Distribuo Normal Modelo",xlab="Valores", ylab= "Quantis da Mostra") # Compara o quartil terico com o quartil da amostra
  qqline(Eco$Temperaturas [Eco$Variables=="Modelo"], col="red")# Mostra a linha de tendncia de normalidade
  qqnorm(Eco$Temperaturas [Eco$Variables=="Sp1"],pch=16, col="black",main="Ajuste da Distribuo Normal Sp1",xlab="Valores", ylab= "Quantis da Mostra") # Compara o quartil terico com o quartil da amostra
  qqline(Eco$Temperaturas [Eco$Variables=="Sp1"], col="red") # Mostra a linha de tendncia de normalidade
  qqnorm(Eco$Temperaturas [Eco$Variables=="Sp2"],pch=16, col="black",main="Ajuste da Distribuo Normal Sp2",xlab="Valores", ylab= "Quantis da Mostra") # Compara o quartil terico com o quartil da amostra
  qqline(Eco$Temperaturas [Eco$Variables=="Sp2"], col="red") # Mostra a linha de tendncia de normalidade
  qqnorm(Eco$Temperaturas [Eco$Variables=="Sp3"],pch=16, col="black",main="Ajuste da Distribuo Normal Sp3",xlab="Valores", ylab= "Quantis da Mostra") # Compara o quartil terico com o quartil da amostra
  qqline(Eco$Temperaturas [Eco$Variables=="Sp3"], col="red") # Mostra a linha de tendncia de normalidade
  x11() # Abre uma janela para grficos 
  par(mfrow=c(1,1))# divide a janela em 1 fila e 1 coluna
  par(bty="l.5", tcl=-0.3) # me permite fazer o fechamento da caixa do grfico, muda a disposio da vinheta nos eixos
  plot (Eco$Temperaturas [Eco$Variables=="Modelo"], pch=16, ylab="Temperatura  C", xlab="Clina de Temperaturas", col="red", ylim=c(20,30))
  boxplot(Eco$Temperaturas ~ Eco$Variables, main="Medidas de tendncia Central", ylab= "Temperatura  C", xlab="Variveis", ylim=c(0,40)) # Grfico de caixas e bigodes
  text (1,35, "a", cex= 1.8) # Inserta a letra "a" na posio 1 no eixo X e 35 no eixo Y
  text (2,35, "a", cex= 1.8) # Inserta a letra "a" na posio 2 no eixo X e 35 no eixo Y
  text (3,35, "a b", cex= 1.8) # Inserta a letra "a e b" na posio 3 no eixo X e 35 no eixo Y
  text (4,35, "b", cex= 1.8) # Inserta a letra "b" na posio 4 no eixo X e 35 no eixo Y
  Variaveis <- tapply (Eco$Temperaturas, INDEX=Eco$Variables, FUN= summary)# Resumo das Variveis
  Sd <- tapply (Eco$Temperaturas, INDEX=Eco$Variables, FUN= sd)# Desvio padro das Variveis
  Var <- tapply (Eco$Temperaturas, INDEX=Eco$Variables, FUN= var)# Varincia das Variveis
  Anova <- lm(Eco$Temperaturas ~ Eco$Variables)# Faz o teste de ANOVA
  termorregulacao <- summary.aov (Anova)# Mostra o resumo da ANOVA
  levene.data <- abs (Anova$residuals) # valor absoluto dos residuos da ANOVA 
  levene <- aov (levene.data ~ Eco$Variables) # Faz o teste de Levene
  termorregulacao2 <- summary (levene) # Mostra o resumo do teste de Levene
  summary<-list ((Variaveis), (Sd),(Var),(termorregulacao), (termorregulacao2) ) # Mostra os resumos de cada variavel, o resumo da ANOVA e o resumo do teste de Levene
  return(summary)
  par(mfrow=c(1,1))
  }
Ecofisiologia (Eco)